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ABSTRACT 

Cen X-4 is one of the best studied low-mass neutron star transients in quiescence. Thanks to XMM- 
Newton large throughput, Cen X-4 was observed at the highest signal to noise ever. This allowed us to 
disclose rapid (> 100 s), large (45 ±7 rms in the 10~^ — 1 Hz range) intensity variability, especially at low 
energies. In order to highlight the cause of this variability, we divided the data into intensity intervals 
and fit the resulting spectra with the canonical model for neutron star transients in quiescence, i.e. an 
absorbed power law plus a neutron star atmosphere. The fit is consistent with a variable column density 
plus variability in (at least) one of the spectral models. Variations in the neutron star atmosphere might 
suggest that accretion onto the neutron star surface is occuring in quiescence; variations in the power law 
tail should support the view of an active millisecond radio pulsar emitting X-rays at the shock between 
a radio pulsar wind and inflowing matter from the companion star. 

Subject headings: accretion, accretion disks — binaries: close — star: individual (Cen X-4) — stars: 
neutron 



tical spectrum shows the characteristics of a K5-7 main 
sequence star, contaminated by lines (e.g. Ha, H/3 and 
H7) and continuum emission probably resulting from an 
accretion disk (Shahbaz, Naylor & Charles 1993). The 
latter was estimated to contribute ^ 80%, ^ 30%, ^ 25% 
and ~ 10% of the quiescent optical flux in the B, V, R 
and I bands, respectively. A short wavelength HST/FOS 
spectrum of the quiescent optical counterpart yielded a 
I35O-22OOA luminosity of 0.6 - 4 x lO^i erg s"^ These 
results were confirmed by higher quality HST/STIS spec- 
tra (McChntock & Remillard 2000). In particular, the UV 
spectrum appears to lie on the extrapolation of the power 
law component seen in X rays. 

In this paper we take advantage of an XMM-Newton ob- 
servation of Cen X-4 to investigate, in much higher detail 
its quiescent state. We describe the observation charac- 
teristics and data filtering in Section 2. In Section 3 we 
discuss on spectral and timing results. The interpretation 
of the data and discussion is presented in Section 4. Our 
conclusions are in Section 5. 

2. OBSERVATION AND DATA ANALYSIS 

The XMM-Newton Observatory (Jansen et al. 2001) 
comprises three ^ 1500 cm^ effective area X-ray telescopes 
each with a European Photon Imaging Camera (EPIC, 
0.1-15 keV) at the focus. Two of the EPIC imaging spec- 
trometers use MOS CCDs (Turner et al. 2001) and one 
uses pn CCDs (Striider et al. 2001). Reflection Grat- 
ing Spectrometers (RGS, 0.35-2.5 keV, den Herder et al. 
2001) are located behind two of the telescopes. 

Cen X-4 was observed on 2001 August 20-21. The EPIC 
cameras were operated in the Prime Full Window mode 
with thin fllters. We used the data generated by the 
Pipeline Processing Subsystem in September 2001. The 
observation is plagued by high background intervals due 
to soft proton flares. We excluded them using an inten- 
sity filter: all events accumulated in a time interval such 

^ INAF-Osservatorio Astronomico di Brera, Via Bianchi 46, 1-23807 Merate (Lc), Italy 
^ INAF-Osservatorio Astronomico di Roma, Via Frascati 33, 1-00040 Monteporzio Catone (Roma), Italy 
CNR-IASF, Istituto di Astrofisica Spaziale e Fisica Cosmica Sezione di Milano "G.Occhialini" , Via Bassini 15, 1-20133 Milano, Italy 
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1. INTRODUCTION 

Many Low Mass X-ray Binaries (LMXRB) accrete mat- 
ter at very high rates, and therefore shine as bright X-ray 
sources, only sporadically. Among these systems are Soft 
X-ray Transients (SXRTs) hosting an old neutron star (for 
a review see Campana et al. 1998a). These systems alter- 
nate periods (weeks to months) of high X-ray luminosity, 
during which they share the same properties of persistent 
LMXRBs, to long (1-few 10 years) intervals of quiescence 
in which the X-ray luminosity drops by up to 5-6 orders 
of magnitude. 

Cen X-4 is one of the best studied SXRTs. Bright X- 
ray outbursts were detected in 1969 and 1979; the source 
remained quiescent since. During the 1979 outburst Cen 
X-4 reached a peak flux Lx ~ 4 x lO'^'^ erg (for a dis- 
tance of d ~ 1.2 kpc, Kaluzienski, Holt & Swank 1980). 
Type I bursts were observed, testifying to the presence of 
an accreting neutron star. Cen X-4 was observed several 
times in quiescence (Asai et al. 1996, 1998; Campana et al. 
1997, 2000; Rutledge et al. 2001). The spectrum was fit 
with a soft thermal component (neutron star atmosphere 
or black body) plus a hard power law with photon index 
in the 1-2 range. Rutledge et al. (2001) reported that the 
0.5-10 keV luminosity decreased by 40 ± 8% in the 5 yr 
between the ASCA and Chandra observations. This vari- 
ability can be attributed to the power-law component; on 
the contrary temperature variations in the thermal com- 
ponent are limited to ^ 10%. Short time variability was 
discovered by Campana et al. (1997) during a ROSAT 
HRI pointing, with the source flux fading by a factor of 
~ 3 on a timescale oi ^ 4 days. 

Extensive spectroscopic and photometric measurements 
of the optical counterpart in quiescence (V=18.7 mag) led 
to the determination of the orbital period (15.1 hr; Cheva- 
her et al. 1989; McClintock & RemiUard 1990; Cowley 
et al. 1988) and mass function (~ 0.2 A/©, converting 
to a neutron star mass between 0.5 — 2.3 Mq). The op- 
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that the rate, excluded the source region, exceeded 0.3 
c in the 10-12 keV band for each of the two MOS 
cameras and 0.8 c s~^ for the pn camera were rejected. 
We adopted these thresholds, which are somewhat higher 
than the standard ones, motivated by the strength of our 
source (~ 1 c s~^ in the 0.2-8 keV band of the pn, see be- 
low). We verified that results obtained by using standard 
thresholds are consistent with the cleaner (but shorter) ex- 
posures. We obtained net (original) exposure times of 32 
(52) ks, 32 (52) ks, 24 (40) ks for MOSl, M0S2 and pn, 
respectively. These represent the deepest observation of 
a neutron star SXRT ever carried out. RGS spectra were 
heavily affected by soft proton flares due to the dispersion 
over the field of view. Optical monitor was not operated. 

For the pn, single and double pixel events were selected 
(patterns to 4), for the MOS, events corresponding to 
patterns from to 12 were selected. Source spectra were 
extracted from a circular region of 45" for the two MOS 
cameras and 40" for pn camera, corresponding to an en- 
circle energy fraction of ~ 90%. Background spectra have 
been extracted from nearby circular regions 2' and 2.5' for 
the MOSl and M0S2, respectively. For the pn we consid- 
ered two circular regions of 40" and 60" in the same CCD 
of the source. Following the prescriptions of Snowden et al. 
(2002) we extracted the spectra (including the FLAG=0 
option) obtaining 9609 counts (0.29 c s~^), 9829 counts 
(0.30 c s-i), 28673 counts (1.1 c s'^) from the MOSl, 
M0S2 and pn in the full band, respectively. Note that 
pile up was negligible (< 1%). 

Response matrices and ancillary region files were gener- 
ated with the SAS (v5.3.3) tasks rmf gen and arf gen. The 
spectral analysis was carried out in the 0.5-10 keV energy 
range for all the instruments using XSPEC (vll.2). All 
spectral uncertainties are given at 90% confidence level 
for one interesting parameter (Ax^ = 2.71). 

3. RESULTS 

3.1. Spectral analysis: I 

We first consider the total spectrum. We rebinned the 

MOS and pn spectra so as to have 25 and 30 counts per 
spectral channel, respectively. We verified that the back- 
ground spectra, especially at high energies, are not domi- 
nated by residual soft proton flares. In particular, pn and 
MOS background spectra show the same behaviour and 
decrease by a factor of ^ 2 from 2 to 10 keV as expected 
for normal conditions (Lumb et al. 2002). 

The spectrum cannot be flt with any single component 
model. We then tried the canonical spectrum for SXRT in 
quiescence (e.g. Campana et al. 1998b; Campana 2001), 
i.e. an absorbed neutron star hydrogen atmosphere plus 
a power law (we adopted the absorption model TBABS by 
Wilms, Allen & McCray 2000 and the neutron star atmo- 
sphere model by Gansicke, Braje & Romani 2002). We 
included also a constant factor to account for the mis- 
match between the different EPIC instruments (Kirsch 
2002). This model provides a reduced Xred ~ ^■^'^ ^'^^ 
518 degrees of freedom (d.o.f.), corresponding to a null hy- 
pothesis probability (n.h.p.) of 12% (see Fig. 1). Neutron 
star atmosphere models with iron or solar composition or a 
disk black body model provided worse flts with Xred ^ l-^- 
A black body model formally provided a better fit to the 
data but with a null column density (< 1.4 x 10^° cm~^ see 



Table 1) . Taking instead a column density of 9 x 10^° cm~^ 
(in line with expectations from optical data) we obtained 
a worse fit (x^ed = 1-18). 

The unabsorbed 0.5-10 keV flux of the power 
law plus neutron star atmosphere models is 2.3 x 
10~^^ erg s~^ cm~^ that at a distance of 1.2 kpc translates 
into a luminosity of 3.9 x 10^^ erg s~^. In the 0.5-10 keV 
band the atmosphere model accounts for 63% of the total 
luminosity. 

3.2. Temporal analysis 

For timing analysis we considered the entire observation, 
including high background time intervals. Besides source 
light curves, we extracted background light curves which 
we subtract (including the error on the background rate) 
from the source light curve. The resulting light curve is 
clearly not constant (Fig. 2). The importance of the back- 
ground is in any case minimal apart from the minimum 
around 21.5 hr in Fig. 2 when the background reached 
40% of the Cen X-4 rate. 

We extracted light curves in selected energy bands: 0.2- 
1 keV, 1-2 keV and 2-8 keV (see Fig. 3). A smaller ampli- 
tude of the variability at high energies is apparent. Taking 
a bin time of 500 s (in order to have a mean of 30 coimts 
per bin in the 2-8 keV light curve) we fit the light curves 
with a constant and we obtained Xred — 

13.7, 4.8, 1.6 

in the three energy bands, respectively. Possible fiare-like 
events can be identified around 14 hr, 16.5 hr and 22.5 hr 
in Fig. 2. Rise times are about 20 min. A dip-like feature 
is also visible around 21.5 hr. 

Wo obtained a power spectrum of the entire observation 
by using the cleanest pn exposure. The power spectrum 
shows a strong noise component below ~ 0.01 Hz (i.e. a 
timcscalc longer than ^ 100 s), consistent with the vari- 
ability discussed above. The rms variability amounts to 
45 ± 7% in the 10~^ — 1 Hz band. The power spectrum 
can be well fit with a power law with index —1.2 ±0.1. No 
significant periodicities or quasi-periodic oscillations are 
seen. The power-law index of the X- ray power spectrum 
is consistent with the optical one which is characterized by 
an index of ^ —1, —1.5 (Hynes et al. 2002; Zurita et al. 
2003). 

3.3. Spectral analysis: II 

To study possible spectral changes related to the vari- 
ability we constructed the color-color diagram made by a 
soft color (1-2 keV / 0.2 1 keV) versus a hard color (2-8 
keV / 1-2 keV, see Fig. 4). Spectral variations are visible 
in the data. We plot in the same diagram the colors ex- 
pected for single (absorbed) power law and neutron star 
atmosphere spectra for a set of parameters. As expected 
from the spectral analysis (section 3.1) these single com- 
ponent models cannot account for the data. In addition, 
we plot two lines with different column densities and dif- 
ferent fractions of the two components models. We note 
that variations appear to be consistent with variations in 
the column density across the observation. Alternatively 
a variation in the fraction of the power law flux to neu- 
tron star atmosphere flux can account for the observed 
variation. 

To check these hypotheses we carry out a more detailed 
spectral analysis. We concentrate on pn data, since these 
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provide a factor of 2 more counts than the two MOS 
cameras together. We divided the data in tlie pn light 
curve (1,000 s bin) at three different count rates: above 1.3 
c s~^, between 1.0 and 1.3 c and below 1.0 c s~^. We 
then fit the three corresponding spectra in order to search 
for possible differences. We consider the same absorbed 
power law plus neutron star atmosphere model. The same 
spectrum for the three count rate spectra is clearly not 
acceptable (x^ed = 2.60). Letting the column density vary 
freely produces an improvement in the fit but not enough 
to make the fit acceptable (Xred ^ ^-5^', the column den- 
sity is smaller at larger count rates). This shows that 
the variability cannot be fully ascribed to a column den- 
sity variation. Better results can be obtained allowing (at 
least) one of the spectral component to vary together with 
the column density^ (Table 2). With an additional param- 
eter we obtained better results, however we cannot decide 
on a statistical basis which of the two fits is better. No 
clear correlations are observable. The column density is 
within the errors for both these new fits. In the case of 
a variable neutron star atmosphere, the temperature in- 
creases with the count rate, in the case of a variable power 
law, the power law photon index steepens as the count rate 
increases. 

4. DISCUSSION 

The quiescent state of Cen X-4 has been recognised to 
be variable both on long times scales (^ 40% in 5 yr, Rut- 
ledge et al. 2001) and on shorter timescales (factor of ^ 3 
in a few days, Campana et al. 1997). During this XMM- 
Newton observation X-ray variability has been observed 
(at a level of ^ 45% rms) on even shorter timescale (down 
to ~ 100 s) thanks to the large collecting area. Variabil- 
ity on such a short timescale would have been missed if 
observed with any previous X-ray satellite. 

In the pn light curve three flare-like events can be iden- 
tified. X ray flares have been observed in the transient 
black hole candidate V404 Cyg also in quiescence (Wag- 
ner et al. 1994; Kong et al. 2002). Flare activity has 
been recently reported also in the optical for a number of 
transient black holes during quiescence as well as for Cen 
X-4 (Zurita, Casares & Shabhaz 2003; Hynes et al. 2002). 
Flares occur on timescales of minutes to a few hours, with 
no dependence on orbital phase. R-band luminosities are 
in the range 10^^ — 10^^ erg s~^. The mean duration of 
optical flares in Cen X-4 is 21 min. This is similar to what 
observed in the X rays. 

Several mechanisms have been proposed for optical 
flares (Zurita et al. 2003; Hynes et al. 2002). Chro- 
mospheric activity, instabilities in the mass transfer rate 
from the companion and viscous instabilities in an accre- 
tion disk have been ruled out for different reasons (Zurita 
et al. 2003). Reprocessing of X-ray variations in the outer 
disk regions has also been considered less likely due to the 



large ratio between optical and X-ray emission (Zurita et 
al. 2003). In passing, we note that in the X ray light curve 
three flare events can be identified (see Fig. 2). The dura- 
tions and energetics of the first two flares are comparable: 
0.9 and 1.2 hr and 1.8 and 4.4 x lO'"*''' erg, respectively. 
The last flare is much shorter and less energetic (0.1 hr 
and 0.6 x 10^^ erg). 

Long term X-ray variability in the quiescent state of 
Aql X-1 has also been reported based on Chandra data. In 
particular, Rutledge et al. (2002) showed that these varia- 
tions cannot be caused by variations in the power law only, 
but can be accounted for by a varying temperature of the 
neutron star atmosphere model. These variations are not 
monotonic but show first a decrease, then an increase and 
finally a decrease in time. This would rule out the deep 
crustal heating mechanism for SXRTs in quiescence. Cam- 
pana & Stella (2003) provided a different interpretation of 
these data (plus an additional BeppoSAX observations) : a 
varying column density and power law component together 
with a stable neutron star atmosphere model. 

This latter model is motivated by the recent discovery 
of a millisecond radio pulsar PSR J1740-5340 in the glob- 
ular cluster NGC 6397 (D'Amico et al. 2001; Ferrario et 
al. 2001). This pulsar shows irregular eclipses in radio and 
it is also emitting in X-rays (Grindlay et al. 2001). The 
X-ray emission is likely powered by the interaction of the 
relativistic wind of the radio pulsar with matter outflow- 
ing from the companion. This shock emission mechanism 
has also been put forward to explain the X-ray emission 
of the radio pulsar PSR 1259-63 orbiting a Be companion 
(Tavani & Arons 1997). This model is not disproved by 
the data. The power law photon index is between 1.5 and 
2 as one would expect from simple modelling based syn- 
chrotron emission. The luminosity increases as the power 
law photon index increases as observed in PSR 1259-63 
for r < 1.8. 

5. CONCLUSIONS 

We observed large (^ 45%) r.m.s. X ray variability in 
the quiescent state of Cen X-4. Small spectral variations 
are observed as well. These can be mainly accounted for 
by a variation in the column density together with another 
spectral parameter. Based on the available spectra we can- 
not prefer a variation of the power law versus a variation in 
the temperature of the atmosphere component (even if the 
flrst is slightly better in terms of reduced x^). This vari- 
ability can be accounted for by accretion onto the neutron 
star surface (e.g. Rutledge et al. 2002) or by the variable 
interaction between the pulsar relativistic wind and mat- 
ter outflowing from the companion in a shock front (e.g. 
Campana fc Stella 2003). 

We thank M. Guainazzi for his advice in the data anal- 
ysis. 



^ In the case of a free neutron star atmosphere component, we consider only a variation in temperature and check that indeed the neutron star 
radius remains stable. 
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Table 1 

Spectral models for the entire XMM-Newton observation. 



Model 


Nh 


kT 


Photon 




n.h.p. 




(lO^o cm-2) 


(eV) 


index 


(d.o.f.) 




Hydrogen 


8.3 ±3.5 


85 ±3 


1.56 ±0.09 


1.07 (518) 


12.0% 


Solar 


11.1 ± 1.5 


> 14 


3.41 ±0.08 


2.88 (518) 


0.0% 


Iron 


4.5 ±1.0 


< 273 


2.79 ±0.09 


2.23 (518) 


0.0% 


Disk BB 


7.3 ±2.2 


219 ± 13 


1.73 ±0.13 


1.31 (518) 


0.0% 


Blackbody 


< 1.4 


186 ±5 


1.64 ±0.09 


1.06 (518) 


16.1% 



Table 2 

Standard spectral models fits for the high and low count rate interval. 



Model 


Nh 


kT 


Photon index 


/(-'red 


N.h.p. 




(1020 cm-2) 


(eV) 


(hard component) 


(d.o.f.) 




Low 


6.0 ±2.9 


87 ± 1 


1.53 ± 0.17 (40%) 


2.60 (450) 


0% 


Medium 








High 












Low 


15.1 ±3.6 


78 ±9 


1.60 ± 0.14 (35%) 


1.51 (448) 


0% 


Medium 


8.9 ±4.3 










High 


4.4 ± 3.4 










Low 


5.1 ±2.6 


86 ±7 


1.47 ±0.13 (46%) 


1.02 (446) 


37% 


Medium 


4.6 ±2.6 


91 ±5 


(38%) 




High 


6.0 ±2.1 


99 ±7 


(30%) 






Low 


9.6 ±4.6 


85 ±6 


1.32 ±0.19 (33%) 


1.03 (444) 


34% 


Medium 


5.4 ±3.7 




1.62 ± 0.16 (42%) 






High 


4.6 ±3.8 




2.11 ±0.20 (53%) 







Note. — The emitting radius of the neutron star atmosphere component has been constrained to be the same in all the three observations. 
In the case of a varying neutron star atmosphere model the radius is 9.7 ± 2.4 km (intrinsic radii), whereas for a varying power law we have a 
radius of 11.5 ± 6.0 km. 
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Fig. 2. — Background subtracted 0.2—8 keV pn light curve of Cen X-4. The bin size is 500 s. Open dots represent the background hght 
curve scaled to the source extraction area (which has been already subtracted in the light curve above). 



7 





> 








1 












K 




0.5 


o 


6 





+ + 



^ 0.3 
cv 0.2 
o ^0.1 

0.3 



> 

^ 0.2 



o 



03 



15 20 

Time (hours) 



25 



Fig. 3. — Background subtracted 0.2—8 kcV pn light curve of Cen X-4 in three energy bands: 0.2—1 keV (upper panel); 1—2 keV (middle 
panel) and 2—8 keV (lower panel). The bin size is 500 s. 
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Fig. 4. — Color-color diagram of the pn background subtracted light curve. The soft color is defined as the ratio of the counts in the 0.2—1 
kcV to 1—2 kcV, the hard color as the ratio of the counts in the 2-8 keV to 1—2 keV. On top of this we depict hardness ratios for an absorbed 
power law model (with photon index from 1.2 to 2) and an absorbed neutron star atmosphere model (with temperatures from 70 to 100 eV) 
for different values of column densities (3, 5, 7, 9 x 10'^'' cm~^, dotted lines). Two different (dashed) curves connect the single component 
models with increasing contribution from the two for column densities of 5 and 9 X 10^" cm~^. 



